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ENERGY SPECTRUM OF THE HYDROGEN ATOM IN A SPACE WITH 
ONE COMPACTIFIED EXTRA DIMENSION, x 


MARTIN BURES 


Abstract. We investigate the consequences of one extra compactified dimension for the energy spectrum of 
the non-relativistic hydrogen atom with a potential defined by Gauss’ law, i.e. proportional to l/|a:p in non- 
compactified 4d space. The calculations were performed numerically by diagonalizing the Hamiltonian in two 
different sets of basis vectors. The energy levels and electron probability density are plotted as a function of the 
compactification radius. The occurrence of several physical effects is discussed and interpreted. 


1. Introduction 


In view of modern theories, aspiring to unify gravity with the remaining forces, the study of higher¬ 
dimensional models becomes very important. Surprisingly enough, in spite of many years of intensive research, 
one of the basic and most interesting questions, namely what happens to ordinary atoms when more than the 
usual three dimensions are considered, seems not fully answered. 

The problem of stability of the hydrogen atom in a space with one additional spatial dimension and a 
potential dehned by Gauss’ law was subject to a recent article [T]. Apart from an extra spatial dimension of 
an infinite extent (resulting topology K^), the authors considered an extra dimension compactified to a circle 
(resulting topology x S”^). The authors proved that the hydrogen atom in a compactified universe is stable for 
compactification radii smaller than the critical value Rent = ao/4, where ag is the Bohr radius. Furthermore, 
it was demonstrated that the system possesses an infinite number of bound states, with an energy extending at 
least to the ground state energy of the hydrogen atom. 

The purpose of the present work is twofold. Firstly, to calculate the approximate energy spectrum as a 
function of the compactiheation radius. Secondly, to discuss various effects that appear in the spectrum and 
in the electron probability density when the compactiheation radius increases. It is particularly important to 
investigate if the calculations of energy levels reveal some indication of divergent behavior of the ground state 
energy at compactiheation radii close to the critical value Rent [I]- 

To the best of our knowledge of existing literature, there is a lack of results in this direction, and this 
important problem is still unresolved. Perhaps the closest attempt was that of [2], where the Bergmann-Frishman 
transformation [3] was used to reformulate the hydrogen problem in terms of a four dimensional harmonic 
oscillator. Then, the calculation of energy corrections was performed, however, only a different, simplihed form 
of the potential than that corresponding to an extra compactihed dimension, was treated eventually. 


2. Hydrogen atom in non-compactified spaces 


Some authors dehne the d-dimensional Coulomb problem through the ^/\x\ potential irrespective of the 
number of extra dimensions, e.g. nisi IS]. Let us emphasize that we deal with the physically more justified 
potential as given by Gauss’ law in d dimensions. Such an approach gives the potential Vd{\x\) ~ and 

the corresponding Schrodinger equation reads 




i’ 




( 1 ) 


where Cd is a d-dimensional charge. This is an utterly different situation. A potential obtained from Gauss’ law, 
when treated in higher dimensional spaces with non-compactified extra dimensions, leads to instability [ZllllSl 
uni nn na IS]- This result remains valid also in the case of four spatial dimensions where, however, a special 
treatment is required. Due to the degree of singularity, the potential energy term l/|a;p can be merged with the 
centrifugal term arising from radial reduction of the Schrodinger equation in spherical coordinates. This results 
in the existence of a critical value for the (dimensionless) coupling parameter Z := 2me\/h?. For Z G [0, Zerit), 
the system is stable without any bound state solutions, for Z > Zerit, the energy of the system is unbounded 
from below Dniis]. Thus, there is no stable hydrogen atom if the extra dimension is infinite. However, if we 
compactify the extra dimension below a certain critical value of the compactification radius, then Gauss’ law 
gives a potential that results in a stable system, as was shown in detail in [1]. 
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3. Hydrogen atom in space with one extra compactified dimension, x 

Following [T] , we consider the underlying space to be x 5"^, where is the fourth dimension X 4 compactified 
in a circle of radius R. The parametrization of by using an angular coordinate 9 via X 4 = R0 then leads 
to the compactified configuration space x (— 7 r, 7 r). The potential of a point charge in four dimensions, 
i.e. —e^/lxp, where 64 is the four-dimensional charge (with unit energy^/^ x length), determines the potential 
in the compactified case: 


V(r,e) = - 


64 


+ R^ {9 — 27rn)2 ’ 


( 2 ) 


where r'^ = x^ + x 2 + x'^ is the radius vector in three dimensions. This expression can be summed and written 
in a closed form as 


V{r,9) = - 


sinh (r / R) 


e 

r cosh (r/ii) — cos 6 * r r — 2 e’"/^ cos 0 -I -1 ’ 

where we choose the relation between three- and four-dimensional charges to be: 

e4/2R = 63 = = T? jma, (4) 

with a being the Bohr radius. This relation comes from the limiting behavior of the potential at large distances. 
For details, we again refer to [T]. Here, we will need the Fourier expansion of the potential (§: 


2e^ 


e’'/-^cos6l- 1 


(3) 


n— — oo 


V{r,9)= y u„(r) =-or V{r,9) = -y cos nd. (5) 

nr* 'Y* nr 


2e2 


TL — l 


Thus, the Hamiltonian of the non-relativistic hydrogen atom in a space with one compactified extra dimension 
X is 


H = H44- 


2mR^d(P 

where is the usual hydrogen atom Hamiltonian, 


-vS' 


-\k\r/R\ke 




iJ3d = -^A3d--. 

2 m r 


( 6 ) 


(7) 


For the following, we shall introduce the units of e^/2a for energy and a for length (so that ff'/2m = 1 and the 
elementary charge is fixed to e} =2). 

4. A TOY MODEL IN COULOMB POTENTIAL PERTURBED BY A YUKAWA POTENTIAL 

Let US first illustrate the computational method on a related example. Because the potential in is given 
by a sum of Yukawa-type terms in the r-variable, we will illustrate the use of the method of Hamiltonian 
diagonalization on the simple case of a Yukawa potential. We will then compare the results to those obtained 
by a direct numerical integration of the radial equation. We consider the Hamiltonian 

H = H34- g-e-^\ ( 8 ) 

r 

with H 34 = —A 3 d — 2/r. Here g > 0 and p, > 0 is the strenght and the range parameter, respectively. 

4.1. Hydrogen atom basis. As a first basis set, we choose the bound states of the hydrogen[^ 

{x\nlm) = Rniir)Yim{9,(l)), I G N, m € {-I,... ,1}, n G {I + l,l + 2,...}, (9) 

where 


— 9 


Y (n+iy. yj "-'-1 y j 


( 10 ) 


are the radial function of the hydrogen atom and Yi^i^) are the spherical harmonics. After some calculations, 
we get for the matrix elements of ® the following expression: 


{n'l'm'\H\nlm) = dii>6, 

where the matrix elements Mn^n';i are given by 

pOO 

( 5 , m ) — 2 ^ / drr‘^RniR 
Jo 


o ^nn' m) 


( 11 ) 




I'n'Z ■ 


/O 

g M 

2 \nn' 


1+2 


r/a 

{n-i-iy.{n' -i-iy 
{n + lyin' + ly. 


/o 


^ — <TX 7 - 2 Z + I 

^n—Z —1 


7-2Z + 1 


(12) 


^We use the following definition of the associated Laguerre polynomials (compatible with Mathematica): = EILo (n—"?) 

Then R^i are orthonormal. 
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with a{fi) = 1/n + 1/n' + /i. To evaluate the integral, we use formula [HI Eqn. (14)]0 


M, 




{n - I - ly.in'- I - l)\ {21 + ly. 


{n + iy.{n'+ ly. a 


2i+2 


n' +1 

i' — I — 1 — k 


-) 

na J 


E 

7i—l — l — k 


1 - 


-) 

na J 


n + I 

n — I—1 — kj 

— l — l — k 


1 - 


(13) 


^k + 2l + 1'' 
k 

Numerical computation of matrix elements via is much faster than using the integral form of |12[ We also 
note that we can study the problem independently for each pair of I, m, because states corresponding to different 


values of these quantum numbers do not mix, as we can see from (111, where Kronecker delta terms in indexes 


Z, I' and m', m' are present. For numerical diagonalization, we truncate the Hamiltonian matrixj^by restricting 
the number of basis vectors 0 - we take n S {^ + 1, Z + 2,..., TV} for some suitable value of N. Diagonalizing 
the truncated matrix gives the approximate ground state energy eigenvalues which we plotted in figure 

4.2. Exponential basis. If we instead use the (non-orthogonal) basis set of exponential functions 

(f|n) = 2a3/2e-“"^ n S {1, 2,..., TV}, (14) 

where {ctn}n=i i® some suitable set of parameters, we get (for I = 0) the matrix elements of the Hamiltonian ([^: 

(2y^ Oi^Oiri) 


{m\H\n) = -{am + «« - aman){m\n) - g 


(c^n 




where (m|n) are the overlap integrals obtained from (14) as 

{m\n) = 4(ama„)^/^ 


drr^" 


{oijn+Ot-n 


2^c 


(15) 


(16) 


We choose the parameter set 

= nel,...,N, (17) 

where we used values A = 0.1, B = 1.5, after optimizing ”by hand”. The ground state energy, obtained as a 
results of Hamiltonian diagonalization, is shown again in Fig. 

4.3. Wronskian method. In this simple case of a Coulomb plus a Yukawa-type interaction in one dimension, 
it is easy to verify the calculations by supplementing them by a method based on the use of Wronskians. Let 
us briefly describe it. The usual substitution R{r) = u{r)/r in the Schrodinger equation corresponding to 
Hamiltonian 0 gives: 

',2 , i(Z + l) 2 2e-'^’" 


u''{r) =[e + 


u{r), 


(18) 


r r 

with the energy eigenvalue E = —k^. For small r, the solution is proportional to the asymptotic behavior 
at large r is c/. ns p.ll8]. This motivates us to define two functions Uin{r),Uout{r), being the solutions 


of 18 and satisfying the initial conditions (considering only I — 0): 

Uout{ri) = ri, M'„(ri) = 1, Uout{r3) = Kutirs) = -kUout{r3)- (19) 

Let us note that the error of the method will be proportional to the error we make by imposing the initial 


0.0 


0.1 


0.2 


0.3 


0.4 


conditions (19) at Hnite distances (and we can estimate the error by varying the initial conditions). 

To hnd the energy E = —k"^ (for a given g) of the ground 
state, we continue as follows: we define the Wronskian 
W{r2) = Uout{r2)u[^{r2) - uLt(^ 2 )um(T 2 ), where u^n{r2), 

Uout{r2) are function values at some intermediate point r 2 
(ri < r2 < Ta), obtained by numerical integration of the 
respective initial value problems. We ”sew” the two solu¬ 
tions Uin{r), Uoutif) by requiring them to be linearly de¬ 
pendent through the condition W{uin{r2),Uout{‘f’2)) = 0. 

The requirement of vanishing of the Wronskian at only 
one point r2 is sufficient, because it implies its vanish¬ 
ing for all r. This can been seen from the fact that 


0.5 


H 


= 0, as can be readily 


and from (18). 



bF {ni^^iiout) — "a 
checked by inserting u'l 

We can see from figure that for larger values of l//i, 
the expansion on hydrogen atom basis does not follow the 
Wronskian method calculation as good as the expansion on 

the exponential functions. This may be an indication that the basis set of exponential functions will be a better 
choice then the hydrogen atom basis also for the studied potential of 0, consisting of a sum of Yukawa terms. 


Figure 1. Toy model: ground state energy as a 
function of 1 /g. 


^ The referred formula is slightly more general. For our purpose, we can simplify it using 1 F 2 (a, fe; 6; z) = (1 — z) “to get 
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5. Energy levels of the hydrogen atom in x 


We shall now use the two basis sets of Secs |4.1| and |4.2| to construct suitable bases for diagonalizing the 
Hamiltonian defined on space x The periodicity of space along the compactified dimensions means 
that the wave function is periodic in the fourth coordinate 9. This suggests to take plane waves as a basis of S^. 

5.f. Basis constructed from the hydrogen atom eigenstates. We now use the basis setj^to construct a 
basis of space x by adding plane waves along the extra dimension: 


{x\nlmq) = 




I € N, TTi € {—/,..., Ti ^ T f, / + 2,... q ^ 'Zi. 


( 20 ) 


Here q is a new quantum number that specifies the particular Fourier mode along the compactified dimension. 
Calculation of the matrix elements from Q and (20) gives: 

..2 


{n'I'm' q'\H\rLlmq) = SuiS^ 


^nn' ^qq' 


1 




-T 2 + ^ ) - (1 - V) \q - q'\/R) 


( 21 ) 


where M„_„/;i(l, \q — q'|/i?) is given by (13) and consequently a{\q — q'\/R) = 1/n + 1/n' + |q — q'\/R- Here 
(1 — Sqq') indicates that the corresponding term is present only for q ^ q'. As in Sec. |4.1[ we can study the 
problem independently for each pair of quantum numbers l,m. For numerical diagonalization, we truncate the 
Hamiltonian matrix by restricting quantum numbers n and q to ranges: 

n G {I + 1,1 + 2,..., N}, q G {—Q, ..., 0,... Q}, (22) 

where N and Q define the number of basis vectors in the truncated basis. We note that, as known from the 
properties of Fourier series, in the limit Q —)■ oo the plane-wave basis becomes a complete basis set (here for 
functions on S^). This is not true for the basis of hydrogen atom bound states which, in the limit N ^ oo, 
remains an incomplete basis set for functions on space (the eigenstates corresponding to the continuous 
spectrum are missing). 

Diagonalizing the truncated Hamiltonian matrix ( |21[ ) numerically using Mathematica, we obtained the ap¬ 
proximate energy eigenvalues which we plotted in figure The lifting of degeneracy of energy levels n = 2,3,4 
for / = 0 and Z = 1 is shown in figure 

The approximate solution of the problem, i. e. the eigenvector obtained as a result of Hamiltonian diagonal¬ 
ization in the (truncated) basis constructed from hydrogen atom eigenstates, is 

^ giijS 

1 plm{x) = 'y' 'V] an,qRnl{r)Yim{^)^=, (23) 

n=lq=-Q 


where an^q are the components of the (normalized) eigenvector (weights of corresponding basis states) and N, Q 
define the size of the truncated basis. The probability density \'tljim{r,9)\'^ of finding the electron at point {r,9) 
is plotted in figure]^ for several compactification radii. The results will be discussed in the concluding Sec. 


5.2. Basis constructed from exponential functions. Let us now modify the basis ( |14[ ). Considering only 
spherically symmetric states (/ = 0 ),let us repeat the calculations by using a (non-orthogonal) set of basis states 
constructed from exponential functions: 

(f|zq) = iG{ 1,2,...,I}, qG{-Q,...,0,...Q} (24) 

V ZTT 

where the parameter set {ai}i is kept as in 0 - The matrix elements of the Hamiltonian in this basis are: 


{jp\H\iq) = 


(jpI*'?) 




(2y4aiQq)'' 


where 


{jp\iq) = 4 (a™a„) 3/2 / drr 2 e-(“-+“^> 


J{q-p)e _ 


q-p\/R)^^ 


27r 


Y Qfi + Otj 


^p,q- 


(25) 


(26) 


are the overlap integrals obtained from ( [24| ). The results of Hamiltonian diagonalization are plotted in figure 
Fig. I^then shows the probability density, obtained as \il>im{r,9)\'^ from an analogical expression to ( [23| ). 

6. Comparison of methods, interpretation, summary 

Let us now analyze the results. In Fig. we plotted the calculated spectrum, as a function of the com¬ 
pactification radius, for the lowest energy levels. The plots on the left-hand side were obtained in basis (20) 


constructed from the bound states functions of the hydrogen atom, whereas the right-hand side refers to ba¬ 


sis (24), constructed from exponential functions. Comparing Figures (A) to (B), or (C) to (D), we can see that 


the exponential basis provides lower energies than the hydrogen atom basis. This confirms our expectations 
from the toy model of a single Yukawa potential, cf. Sec. There, the calculations in the hydrogen atom basis 
could not follow the curve obtained from the Wronskian method as good as the exponential basis. 

Energy levels with quantum numbers q = 0 correspond to electron states which do not possess any 9- 
dependence of the wave function (their wave functions are constant in the compactified direction). Levels with 
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q> 1 are the excited Kaluza-Klein states. As the compactification radius R increases, their energy decreases 
as 1 /R^ 


cf. ( |21| ) or (25). For R large enough, their energy falls down to a range comparable to the bound 
spectrum and they cross levels with g = 0, cf. plots (A-D) in Fig. We note that on the right-hand side of 
plots (B) and (D), some of the levels seem to change direction abruptly. This is a mere consequence of the fact 
that for each R we only plot a certain (fixed) number of lowest eigenvalues, hence the continuations after the 
’’break” are not the same levels but the excited Kaluza-Klein states (with g > 1, see above). 

It was shown in [1] that the hydrogen atom with an additional compactified dimension becomes unstable 
for Rcrit > a/4. Calculations using basis (20) give for R = a/4 the ground state energy Ei = —1.00665, the 
exponential basis gives Ei = —1.08315 (both results for = 10, Q = 30). Figs.[^ parts (E) and (F), show the 
dependence of the ground state energy on the respective bases sizes. When truncating to a finite basis, we do 
not get a clear sign of divergence for compactification radii close to the critical value Rcrit- However, we have 
confirmed the that the spectrum changes dramatically, when we approach the critical compactification radius 
Rcrit- For small compactification radii, we can trust our results which agree between the different methods and 
do not depend on the size of the basis. For a more precise calculation of the spectrum close to Rcrit one would 
have to use a different technique. 

In Fig. we can see how the presence of an extra dimensions lifts the degeneracy of levels n = 2, 3 for quantum 
numbers I — 0 and I = 1 (the n = 1, g = 1 Kaluza-Klein state is also visible). This is an effect analogous to 
the splitting of the spectral lines of atoms due to various corrections (e.g. fine or hyperfine structure). Recent 
progress in spectroscopic methods could serve as a testing ground providing further constraints on the size of 
extra dimensions HZl, as the effects predicted by theory can be directly measured: comparing the calculated 
spectrum to experimental values would restrict the possible values of the compactification radius R. 

Finally, Fig.j^shows the electron probability density, as calculated by using the two sets of basis vectors. The 
top row corresponds to basis (20) constructed from bound states functions of the hydrogen atom, whereas the 
bottom row is based on calculations in basis (24), constructed from exponential functions. We can observe how 
the electron tends to localize at the nucleus, as the compactification radius increases. To achieve localization, 
Fourier modes along the compactified direction are required. However these terms increase the energy of the 
electron (which is proportional to 1/i?^, cf. (211 or (25)), so for a given value of R, some sort of balance has to 
be found. Again, for large enough compactification radius, close to Rcrit = a/4, cf. [1], complete localization or 
the eventual ’’fall to the center”, cf. [m p.53-54,114-177], is expected. 

The present calculational method (partial exact diagonalization) provides an upper bound on the real energy. 
Hence, it would be of great interest to use some other method to obtain also a lower bound on the spectrum, 
and to better explore the behavior of the spectrum close to the critical compactification radius. 
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(a) Energy levels: hydrogen atom basis 
{I = m = 0, size N = 10, Q = 30) 


(b) Energy levels: exponential basis 
(Z = m = 0, size A'’ = 10, Q = 30) 
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(c) Energy levels: hydrogen atom basis 
{I = m = 0, size N = 10, Q = 30) 


(d) Energy levels: exponential basis 
{I = m = 0, size A'" = 10, Q = 30) 
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(e) Ground state energy dependence on basis size 
(hydrogen atom basis, N = 7, Q = 1,..., 50) 
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(f) Ground state energy dependence on basis size 
(exponential basis, N = 7, Q = 1,..., 50) 


Figure 2. Energy eigenvalues (in units j2a) as a function of the compactification radius R 
(in units of the Bohr radius a). Hydrogen atom basis (left-hand side), exponential basis (right- 
hand side). The computational step in R was adjusted according to the second derivative of 
the curves between AR = 0.005 and AR = 0.03. 
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Figure 3. Lifting of degeneracy of energy levels (hydrogen atom basis, N = 7, Q = 30): 
n = {2,3}: I = 0 (solid line) I = 1 (dashed line), m = 0. The almost vertical curve represents 
the hrst Kaluza-Klein state n = l,q = 1. 
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(a) Hydrogen atom basis {N — 10, Q = 30) 



(b) Exponential basis {N = 10, Q = 30) 

Figure 4. Probability density in the {r,6) plane for several values of R. 





















































